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In this work we study the spectral functions of scalar mesons in one- and two-channel cases. 
When the propagators satisfy the Kallen-Lehman representation a normalized spectral function is 
obtained, allowing to take into account finite-width effects in the evaluation of decay rates. In 
the one-channel case, suitable to the light sigma and k mesons, the spectral function can deviate 
consistently from a Breit-Wigner shape. In the two-channel case with one subthreshold channel the 
evaluated spectral function is well approximated by a Flatte' distribution; when applying the study 
to the ao (980) and /o(980) mesons the tree- level forbidden KK decay is analysed. 



INTRODUCTION 



The scalar mesons below 2 GeV are in center of debate since many years [l], 0, 5L 0] • More states than expected 
from the quark-antiquark assignment are reported in Particle Data Group (PDG) 5], leading to the introduction of 
a scalar glueball 0, tetraquark states Q and mesonic molecules In particular, the scalar resonances below 1 
GeV have appealing characteristics, such as the reversed level ordering of masses, expected from tetraquark states 
, Oil QjJ • In turn, this scenario implies that quarkonia lie between 1 and 2 GeV. A complication in the analysis 
of scalar states is mixing: between 1-2 GeV a quarkonia-glueball mixing in the isoscalar sector is considered, for 
instance,in Refs. [Hj]. Mixing among tetraquark states below 1 GeV and quarkonia above 2 GeV is studied in Refs. 
[HI, 0, EH, where, however, the results do not coincide: while a large mixing is found in Ref. [l4j, a negligible 
JL | mixing is the outcome of [l5j]. It should be stressed that different interpretations of scalar state are possible: a nonet 
qj- of scalar quarkonia is settled below 1 GeV in Refs. [l6[ in agreement with the linear sigma model and the Nambu 
Jona-Lasinio (NJL) model, while in Ref. [13] a broad glueball, to be identified with /o(600), is proposed. We refer to 
1—1 . Refs. [2, 0, & E3] for a discussion of arguments in favour and/or against the outlined assignments. 
_ -. ' Studies on scalar mesons have been extensively performed by using chiral perturbation theory [18j, where a scalar 
resonance at about 440 MeV is inferred out of pion-pion scattering. A full nonet of molecular-like scalar states is 
■^j- ■ generated in the unitarized chiral perturbation theory of Ref. (l9| . In particular, Pelaez [2(| studied the large- N c 
, dependence of the light scalar resonances finding that they do not scale as q uarkonia but in agreement with a molecular 
l/") ■ or tetraquark composition (see, however, also the discussions in Refs. [2ll| ) . 

In the present paper we concentrate on an important aspect of light scalar resonances, namely the form of their spec- 
. tral functions, in a simple theoretical context. In this study, relevant for both quarkonium or tetraquark assignment 
of light scalars [22j , effects arising from loops of pseudoscalar mesons are considered: this leads to parametrizations of 
spectral functions beyond the (usually employed) Breit-Wigner and Flatte distributions and allow to include finite- 
width effects in the evaluation of decay rates. In particular, we consider the following physical scenarios: (i) the 
case of a broad scalar resonance, strongly coupled to one decay channel, such as the a = /o(600) in the pion-pion 
decay mode, for which the spectral function can deviate substantially form the Breit-Wigner form; (ii) the case of 
two channels, one of which is sub-threshold for the mass and thus forbidden at tree-level, as the KK decay mode for 
the resonances /o(980) and ao(980). In the latter case a comparison with the usually employed Flatte distribution is 
performed. 

The key-quantity of the discussion is the propagator of scalar resonances dressed by mesonic loops in one or more 
channels. When the Kallen-Lehmann representation is satisfied, as verified at one-loop level in the case of light 
scalar mesons for large ranges of parameters 22], the spectral function (proportional to the imaginary part of the 
propagator) is correctly normalized and is interpreted as a 'mass distribution' for the scalar state. A general definition 
of the decay of a scalar state, involving the obtained mass distribution, is then possible. In this way one takes into 
account in a consistent fashion finite-width effects for the decay, hence allowing to study deviations form the usually 
employed tree-level formula for decay rates. Furthermore, the fulfillment of the Kallen-Lehmann representation offers 
a criterion to delimit the validity of our one-loop study: as soon as violations appear (generally for large coupling 
constants) the obtained spectral functions are no longer usable. 

We regularize the mesonic loops by means of a cutoff function which in turn is equivalent to a nonlocal interaction 
Lagrangian. In a phenomenological perspective it is reasonable that the mesonic states in the loop cannot have 
indefinitely high virtual momenta which are naturally limited due to the finite range of the meson-meson interaction. 
We also show that the dependence on the chosen cutoff function and on the specific value of the cutoff is mild. 

In order to render the paper easily understandable and self-contained we start in Section II with one-channel case 
by recalling the basic definitions and properties, then we apply the study to the scalar sigma a and kaon k resonances: 
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the corresponding spectral function shows consistent deviations from the usual Breit-Wigner one. In Section III we 
turn to the two-channel case, with particular attention to the resonances ao(980) and /o(980), their decay rates and 
spectral functions in comparison with the Flatte distribution [23|, [24| . Implications of the results in view of a nonet 
of tetraquark states below 1 GeV is discussed. In section IV we drive our conclusions, emphasizing as in Ref. [22j 
that the use of propagators fulfilling the Kallen-Lehmann representation, which implies normalized distributions and 
a correct definition of decay rates, should be preferable both in theoretical and experimental works. 



II. SCALAR SPECTRAL FUNCTION: ONE-CHANNEL CASE 



A. Definitions and properties 



We consider the scalar fields S and ip described by the Lagrangian 



4 = \(d,Sf - \m*s* + \{d^) 2 - imV +flV- (1) 
In the limit g = the propagator of the field S reads 

As M = pZ-M> + is > 9 = - (2) 

We intend to study the modification to A${p) when g ^ 0, which arises by considering the loop-diagram of Fig. 1 
and how this contribution affects the decay mechanism S — ► yxp. We recall that at tree-level the decay width reads 



?sL(M ) = ^^[g Svv ] 2 9(M - 2m) , p Sw = J ^° - m 2 , g Sw = V2g. (3) 



8ttM, 







A 2 



where 0(x) is the step function and 

(The factor \[2 in the amplitude #s^ w takes into account that the final state consists of two identical particles.) 
In general, under the symbol psab the expression 



Psab = 2^ + < M A ^ A/ s) 2 ~ 2(Afi + M 2 )M 2 (4) 

i.e. the momentum of the outgoing particle(s), is understood. 

At tree-level the particle S is treated as stable. However, the very fact that the decay TgL ^ for Mo > 2m 
means that S is not stable and cannot be considered as an asymptotic state of the Lagrangian The tree-level 
expression TgL is only valid in the limit g — > 0. The evaluation of the loop of Fig. 1 offers a way to define and 
interpret the decay S — > (f(p as we describe in the following. The modified propagator of S is obtained by (re)summing 
the loop-diagrams of Fig. 1: 



P 2 -M$+g 2 s n P 2 ) + ie 



where the self-energy T,(jr) reads 



^ P } l j (2tt) 4 [(q+p/2) 2 -m 2 +ie}[(q-p/2) 2 -m 2 + ie}' ^ 

The integral defining S(p 2 ) is, as known, logarithmic divergent. Our intention is to consider the Lagrangian £g as an 
effective low-energy description of the fields S and ip, and not as a fundamental theory valid up to indefinitely high 
mass scales. We do not apply the renormalization scheme to but we introduce a regularization function f&(q) 
which depends on a cut-off scale A for the large momenta. The self-energy S(p 2 ) is then modified to: 

Efo* = -i ( -^L jl(gVg) (7) 

U ' J (2tt) 4 [(q + p/2) 2 -m 2 + ie] [(q - p/2) 2 - m 2 + ie] 1 ' 

When choosing /a(<z) = f\(q 2 ) the covariance of the theory is preserved, otherwise is lost. Indeed, in many calculations 
related to mesonic loops a regularization of the kind f\(q) — /a(~9* ) is chosen, which leads to simple expressions 
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FIG. 1: Mesonic loop 



for the self-energy contribution but breaks covariance explicitly, and thus is strictly valid only in the rest frame of 
the decaying particle. In particular, the 3d-cutoff fty{~c? 2 ) = #(A 2 — if 2 ) is often used. We refer to Appendix A for 
a closer analysis of the self-energy S(p 2 ), where the case of unequal masses circulating in the loop is also presented. 
The interaction strength among light mesons is suppressed for distances larger than I ~ 0.5-1 fm: in this particular 
physical example it is then natural to implement a cutoff A ~ l/l, which varies between 1 and 2 GeV . 

The cutoff- function /a(<z) is not present in the Lagrangian Cg of Eq. (fTJ). In this sense the Lagrangian is incomplete 
because it does not specify how to cut the high momenta. One can take into account /a(<?) already at the Lagrangian 
level by rendering the interaction term nonlocal: 

(4) mt = g V - gS J dV(z + y/2M* - y/2)*(y). (8) 

The Feynman rule for the 3-leg vertex is modified as: 

i9 -> ^/A(^y^), / A (ff) = / d 4 yHv)e- ivq (9) 

where q\ and qi are the momenta of the two particles <p. The function /a(<?) enters directly into the expression of all 
amplitudes. In particular, the self-energy contribution of Eq. ([7]) is now obtained by application of the (modified) 
Feynman rules to the loop diagram of Fig. 1. Indeed, the derealization of the interaction term also induces a change 
of the tree-level result for the decay, which becomes (for the case / A (g) = /a("9* 2 )): 

it^(Mo) = |^[35 W /a(Y 2 = P%r)MMo - 2m), g Slpip = ^2g (10) 

that is the function /a("5* 2 ) is explicitly present in the tree-level decay expression and can be interpreted as a 
phenomenological form factor [271 ] 

If a step-function is used the local tree-level expression of Eq. © is recovered, provided that the cutoff A is large 
enough. In this work we use the following cutoff-function 



fA(q) = fA(lt 2 ) = {l + t 2 /^Y 1 - (11) 

With this choice the Fourier transform $(y), see Eq. ©, takes the form 6(y°) exp [— \lf | A] / |"^| , thus decreasing 
rapidly for increasing distance of the two interacting mesons ip. The interaction range I is of the order A -1 , as discussed 
above based on general dimensional grounds. At each step of the forthcoming study we employed also different forms 
of /a(<?); finding that the dependence on the precise form of f\(q) affects only slightly the results. Notice that in Refs. 
[25l ] a similar equation to (fit)]) (where s = if 2 and /a = G(s) in the notation of [25[) represents the starting point of 
the analysis. The function G(s) in the above cited works is taken to be a Gaussian, A is of the order of 1 GeV. The 
present approach shows the link between such form factor /a = G(s) and a nonlocal Lagrangian. However, we will 
not concentrate as in (25| on scattering amplitudes but on spectral functions and decay widths. At the same time we 
do not relate the imaginary and real part of the propagator via Kallen-Lehmann dispersion relation, but we evaluate 
them independently and subsequently we check numerically if it satisfied, see details in the following discussion. 
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FIG. 2: Real and imaginary part of the mesonic loop for m = 0.5 GeV and A = 1.5 GeV. 



Let us now turn to the self-energy £(p 2 ). A general property for S(p 2 ) follows from the optical theorem: 

Hx) = 9s^Irn[E(x = yffi)] = xVf^{x). (12) 

The imaginary part of the self-energy diagram is zero for < x < 2m and nonzero starting at threshold. The real 
part 



R(x)=g 2 Svv Re[E(x=V^)] 



(13) 



is nonzero below and above threshold. In Fig. 2 the functions R(x) and I(x) are plotted using Eq. (fTTj) . A particular 
choice for the parameters m — 0.5 GeV and A = 1.5 GeV (of the order of physical cases studied later) is done. Anyway 
the plotted functions are qualitatively similar for large ranges of parameters. As noticeable, R(x) is continuous but 
not derivable in x — 2m. It has a cusp at x = 2m: the left derivative is +00, while the right derivative is finite and 
negative. 

In terms of the two functions R(x) and I(x) the propagator of Eq. ((2|) reads 



As(x) 



1 



x 2 - M 2 + R(x) + il(x) + is ' 
We define the (Breit-Wigner) mass M for the scalar field S as the solution of the equation 

M 2 - M 2 + R{M) = 0. 



(14) 



(15) 



When the function R(M) is positive, which is usually the physical case (Fig. 2), the dressed mass M is smaller than 
the bare mass Mq, showing that the quantum fluctuations tend to lower it. 

We now turn to the spectral function ds (x) of the scalar field S related to the imaginary part of the propagator as 



2x 

ds(x) — — lim Im[As(x)] 



(16) 



In the limit g — > we obtain the desired spectral function d$(x) — 8(x — Mq). The normalization of ds{x) holds for 
each g 



ds{x)dx = 1. 



The latter Eq. is a consequence of the Kallen-Lchmann representation 

A s (x) = 



2y -Im[A s {y)} 

d v— -j— " 



7T y 



is 



(17) 



(18) 




FIG. 3: Spectral functions in the cases M < 2m (dashed line) and M > 2m (continuous line). The coupling constant is 
gsifrip = 3 GeV and two values of the masses are chosen: M = 0.9 GeV and M = 1.3 GeV, m = 0.5 GeV. 



when taking the limit x — > oo. Eqs. ([T7]) and (fT5|) hold in general for the full propagator. In our case we check 
numerically the validity of the normalization condition (fT7|) at one-loop level of Fig. 1 . We find that it is fulfilled to 
a high level of accuracy for large ranges of parameters, see also the discussion in [22| and in the next subsection. 
Let us consider ds{x) in the two interesting cases M < 2m and M > 2m. If M < 2m Eq. (JTH]) becomes: 

ds(x)=Z6(x-M)6(2m-x) + — — 9 (19) 

7T (x 2 - M 2 + R(x)f + I{x) 2 



where 



2M V dx 



x=M 



-1 



^=(l + 4f¥) ) • (20) 



When M < 2m the constant Z is usually reabsorbed into the definition of the wave function renormalization, hence 
recovering the free propagator properly normalized as (p 2 — M 2 +iej , corresponding to d$(x) = S(x — M) for 
x < 2m as in the free case g — 0. Thus, we have still a stable particle with dressed mass M instead of Mq. Notice 
that < Z < 1 because R'(M) is a positive number: the quantity (1 — Z) can be interpreted as the amount of virtual 
clouds of 2ip contributing to the wave-function. 
If M > 2m the spectral function reads 

, , n 2x I(x) 

d s (x) = ^-L— . 21 

t {x 2 ~ M 2 + R(x)f +I{x) 2 

No delta-functions are present but typically a picked distribution dg(x) is obtained, corresponding to a physical 
resonance. The mass M is not the maximum of the ds(x) although in general very close to it. Consistent deviations 
can appear when M is close to threshold and for large coupling constant, see next subsection for a more detailed 
discussion of this point. Notice moreover that ds(x) is zero for x < 2m. 

We plot the typical behavior of the spectral function in both cases M < 2m and M > 2m in Fig. [3l We used the 
values M = 0.9 GeV and M = 1.3 GeV corresponding to the the two cases below and above the threshold, m = 0.5 
GeV as before and gs<p<p = 3 GeV. The value of Z, for the sub-threshold case, is ~ 0.9, we have numerically verified 
that the spectral functions are normalized in both cases. 

When M > 2m the function ds{x) can be interpreted as the mass distribution of the resonance, see also Appendix 
B for an intuitive discussion about this point. We then define the decay rate for the process S — > tp<p by implementing 
the distribution ds(x), and thus including finite width effects, as: 

/>OC 

T Sw = / dxd s (x)Tf w (x). (22) 
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This formula reduces to the tree- level amplitude r]£* (Mo) of Eq. (fTOf in the limit of small g: 

Tf^Mo) ~ T Sipv for g 0. (23) 

Notice that in this limit M — » Mq. However, even for finite g, when A/ ^ Mo, the formula Ts w ~ r^J, (M) offers a 
first approximation to the decay width of the state as long as the distribution is picked, i.e. the scalar state S is not 
too broad. 

The definition Eq. (|22[) for the decay S — * (pip is thus a generalization of the tree-level result of Eq. (|10p and takes 
automatically into account that the state S has a finite width parametrized by the mass distribution ds{x), which 
naturally arises by considering the self-energy of the scalar propagator. Notice that the real part of the propagator is 
necessary in order Eq. (|17|) to hold: its neglection would spoil the correct normalization. 

Evaluating the real and imaginary part at x = M and neglecting their ^-dependence, the distribution (|21|) is 
approximated by 



4^)-- ^ S (24) 

7T (a; 2 - M 2 ) + (I(M)) 

which is the relativistic Breit-Wigner distribution for the resonance S, usually employed in theoretical and experi- 
mental studies. However, the distribution d)g w (x) neglects the real part of the loop diagram and consequently the 
normalization of Eq. (jTTJ) does not hold, implying that d)g w (x) has to be normalized by hand. At the same time the 
mass M does not coincide with the maximum of ds{x). Thus, we insist on that the usage of automatically normalized 
distribution emerging from propagators fulfilling Kallen-Lehmann should be preferable. 



B. Application to the light scalar mesons a and k 

An interesting example for the one-channel case is the decay of the scalar meson a = Jo (600). As reported by the 
PDG Q, experimental data are affected by large uncertainties both for the value of the mass, M a — 0.4-1.2 GeV, 
and the value of the Breit-Wigner width, = 0.6-1 GeV. The dominant channel, which we will consider here, is the 
decay into two pions, for which M > 2m . 

By applying the formulas of Section II. A we show in the left panel of Fig. [5] the spectral functions d a {x) of the 
cr-resonance for the two boundary cases of PDG, namely M a — 0.4 GeV and M a = 1.2 GeV, respectively, for the 
coupling constant g a7TT! = 3 GeV. 

The spectral function assumes different shapes for different values of the mass. While for M a — 1.2, far from the 
threshold, the spectral function has a regular 'Breit-Wigner-like' form, in the case M a = 0.4 GeV a distorted shape, 
with a narrow peak just above threshold, is visible 28]. The employed value of the coupling constant, g^ix = 3 GeV, 
serves as illustration and actually corresponds to a somewhat too narrow width. The increase of g a7r7r leads, however, 
outside the range of validity of the normalization of Eq. (fTTl) at one-loop level, see below. 

The description of the scalar kaonic resonance k follows the same line [26[ . As shown in the right panel of Fig. @] 
a strong deviation from the Breit-Wigner form is obtained when g^K = 6 GeV, corresponding to Tt = 0.38 GeV, 
while a less distorted shape is found for g^K = 3 GeV, for which =0.11 GeV. 

At this point a short discussion on the definition of the mass of an unstable particle is needed: in the Breit-Wigner 
scheme, the mass of the particle is the value corresponding to the maximum of the spectral function and it is one of 
the parameter of the distribution (the second one is of course the width). In our scheme, using the spectral functions 
coming from the loop evaluation this is no more the case: the mass M defined in Eq. (fT5)) is again a parameter of the 
distribution (together with the coupling constant g) but it does not coincide with the maximum of the distribution. 
While for M a = 1.2 GeV (far from threshold) the maximum of the spectral function occurs at 1.202 GeV, thus only 
slightly shifted from the mass, when M a = 0.4 GeV the maximum of the spectral function (apart form the threshold 
enhancement peak) occurs at sizeable larger values respect to the mass, here 0.454 GeV(29|. (Notice also that in the 
latter case the bare mass Mq is 0.614 GeV, thus implying a strong influence of the pion loop to the sigma mass, see 
also Appendix A for a comparison of different 'masses'). Indeed, although the mass M being the zero of the real 
part of the inverse propagator, see Eq. (fT5")) . is referred to as Breit-Wigner mass 25], the best fit to the full spectral 
function d a (x) by using a Breit-Wigner form is obtained for a Breit-Wigner mass Mbw coinciding with the maximum 
of the distribution. It is also remarkable that in some cases the spectral function has not a maximum, a part from 
the threshold enhancement peak, as we can see for the k meson in the right panel of Fig. [4] for g^K = 6 GeV. 

We now study closer the decay process a = /o(600) — > tttt using Eq. (f2"2")l which implements the spectral function 
d a (x). In Fig. Owe compare the full and the tree-level decay rates for different values of the cut-off and of the coupling 
constants (a mass M a = 600 MeV is used). 
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x [GeV] 
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FIG. 4: Spectral functions of a and k. 

Left panel: A = 1.5 GeV, g = 3 GeV, M a = 0.4 GeV with a correspondent F = 0.242 GeV and M CT = 1.2 GeV with a 
correspondent F = 0.113 GeV. 

Right panel: Mt = 0.8 GeV, A = 1.5 GeV, g = 3 GeV with a correspondent T — 0.112 GeV and g = 6 GeV with a correspondent 
F = 0.382 GeV 

r [GeV] 




g[GeV] 

1.5 2 2.5 3 3.5 4 



FIG. 5: Fhe full (solid) and the tree-level (dashed) decay rates a — > nn are shown as function of the coupling constant. The 
cases A = 1 and A = 2 GeV correspond to thin and thick lines respectively. 

As expected, there is not a strong dependence on the choice of the cut-off. Moreover the results obtained with our 
formulae are well in agreement with the tree-level results for small values of g (since the spectral function tends to a 
delta function). Non- negligible differences instead occur for the larger values of g. This is due to a different analytic 
dependence of the two corresponding formulae from the coupling constant: while the tree-level expression depends 
quadratically on g, in the loop formula g appears also in the distribution ds(x). 

The limit of the validity of the employed one-loop level analysis is an important aspect which deserves further 
discussion. Namely, when the coupling constant is too large and the mass is not far from threshold the normalization 
condition of Eq. (flTf is lost, see also the corresponding discussion in [22j]. This fact means that higher orders must be 
taken into account to satisfy the Kallen-Lehmann representation and thus to recover the correct normalization of Eq. 
(|17[) . At the same time the violation of the normalization is a valid criterion to establish the limit of our study: for 
this reason in Fig. 5 we stop the plot at g a ^n — 4 GeV (corresponding to Y a ~ 400 MeV) , in fact larger values imply 
J °° ds{x)dx < 1. At this point the full decay width (using Eq. ([22]) ) is already ~_60 MeV smaller than the tree-level 
counterpart. A decay width of about 400 MeV is on the low side for the a (see [U)- In Ref. [H} a width 150 MeV 
larger is obtained. A study beyond one-loop level would then be necessary to evaluate the spectral function for larger 
coupling (i.e. larger width) and represents a possible outlook of the present work. Surely the overestimation of the 
tree-level formula keeps growing for increasing interaction strengths. Similar considerations hold for the k meson. 
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III. SCALAR SPECTRAL FUNCTION: TWO-CHANNEL CASE 



A. Definitions and properties 

We now consider two channels for the scalar resonance S described by the Lagrangian density (m 2 > mi): 

C% = \(d,S) 2 - 1 -M 2 S 2 + 1 -(d^ 1 ) 2 - 1 -rn 2 M + 

-^{d^f ~ ^mjipj +giS(fl + g 2 Sifl. (25) 
The processes S — > ifi^pi and S — ► tp 2 (p 2 correspond to the tree- level decay rates 

TsU^Mo) = ^$-[9s<p lV1 ?0(M - 2m!), g SviVl = V2g 1: (26) 
rf V2V2 (Mo) = f^p[2s^] 2 W - 2m 2 ), 9s ^ 2 = V2g 2 . (27) 



'0 



The propagator is modified by loops of (p\ and <p 2 , denoted as £i(p 2 ) and T, 2 (p 2 ) and given by Eq. ([7]) for m = mi 
and m = TO2 respectively. A delocalization of the interaction, via a vertex- function <&(?/) and the corresponding 
Fourier-transform f\(q) = J d i y&(y)e~ lyq , is then introduced as in Eq. |[8j for both channels in order to regularize 
the self-energy contributions. As a consequence, the tree-level results are modified as 

IsU^Mo) = U^[5s^ 1 ^ 1 /A(^ 2 =p| VlVl )] 2 e(Afo-2m 1 ) (28) 



n 



and the propagator as 



IsU 2 (M ) = §^[3Sv, 2 ^/a(^ 2 = p 2 s^ 2 )] 2 0(M - 2m 2 ), (29) 



As(2;) ~ x 2 ~ Af 2 + i?(x) + il(x) + is ■ (30) 



where 



and 



Ri?) = ^^[E^x = VP 2 )] + g 2 Sv2V2 Re[Z 2 (x = vV)] (31) 

m = g 2 SviVl I™[^i(x = ^/7)}+9 2 Stp2lp2 Im[X 2 (x = y/p)] (32) 
= xTf^^+xTf^ix). (33) 

In the last equation the optical theorem has been used. The mass M of the state 5 is given by M 2 — Mg + R(M) = 0. 
Again, we have two cases: 

(i) M < 2m\. the distribution ds(x) takes the form Z5(x — M) for x < 2m\. The discussion is similar to the 
one-channel case. Eq. (fl"9|) is still valid. At threshold 2m\ the continuum starts. 

(ii) M > 2mi: as in Eq. (|21"j) the distribution is 

/ x 2a; I(x) , . 

d s (a;) = . 34 

t (x 2 - M 2 + i?(x)) 2 + /(a;) 2 

It vanishes for M < 2mi. At x = 2m 2 the second channel opens. 

In the case (ii) we have a resonant state. The decay rates into the two channels S — > (fiipi and S — > ^2^2 are given 
by the integrals: 

Ts^^! = / dxd s {x)Y t s l VlVl {x), T SV21P2 = / tfod s (x)r^ 2!p3 (a;). (35) 
Jo Jo 

A particularly interesting case takes place when 2mi < M < 2m 2 : while the tree-level result for S — > (f2f2 vanishes, 
we find that Tstp 2V > 2 is not zero. In this case the tree-level approximation is absolutely not applicable: the particle 
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FIG. 6: Left panel: spectral function of fo within our formalism (solid line) and using the Flatte' distribution (dashed line). 
Parameters are: A = 1.5 GeV, 9f ()K j^ = 3 GeV with a correspondent total decay rate F = 0.058 GeV. 

Right panel: spectral function of ao within our formalism (solid line) and using the Flatte' distribution (dashed line). Parameters 
are: A = 1.5 GeV, g f oK j^ = 3 GeV, with a correspondent total decay rate F = 0.048 GeV. 

S does decay in virtue of the high-mass tail of its distribution. A physical example is well-known: the resonances 
/o(980) and ao(980) have a non-zero decay rate into KK, although their masses are below the threshold Ittik- Clearly, 
a sizable decay rate Ts V2 <p2 1S obtained only when M is close to threshold. A generalization to the present definitions 
to N channels is straightforward [22| . 

When applying the decay formulas (|35p it is however important to verify numerically that the normalization of the 
distribution ds(x) holds: in fact, as discussed in Section II. B, only in this case the formalism is self-consistent. 



B. Application to a (980) and / (980) 

In this subsection we study the spectral functions of the scalar mesons a = a (980) and / = / (980), whose 
masses are M ao = 984. 7± 1.2 MeV and Mf = 980± 10 MeV [f|. For both resonances two decays have been observed: 
ao — > ttt], ao — > KK and fo — > 7T7t, /q — > KK. Notice that both masses are below threshold of kaon-antikaon 
production, M ao j < 2Mjc — 987.3 MeV, thus the decay of both resonances in KK vanishes at tree- level, while 
experimentally it was seen for both ao and /o states. 

For definiteness we use the following ratios obtained in the experimental analysis of [J| 



-!p± = 4.21 ± 0.46, = 2.15 ±0.40, = 0.75 ± 0.11 (36) 

^/o 7r7r ®a KK ^a Q KK 

therefore leaving us with only one free parameter, chosen to be g f oK -^ ■ Although experimental uncertainties are still 
large, the results of Eq. (|3"6")) are qualitative similar to various studies, see [24| and Refs. therein, pointing to a large 
KK coupling for both resonances with a particular enhancement for /o (see 0, OH EH EH and Refs. therein for 
spectroscopic interpretations) . 

For the typical value g f QK -j^ = 3 GeV we report in Fig. [5] the spectral functions of a and f There is a large 
probability, ~ 50%, in both cases, that these two mesons have a mass larger than the threshold of production 
2rriK and therefore the tree-level forbidden decay occurs. In the same figure we compare our distribution with the 
Flatte one (23l . [24| . which is usually employed for the ao and /o mesons. At variance from our distribution, Flatte 
distribution must be normalized by hand. The two distributions are quite similar, only for the fo the values of the 
mass corresponding to the maximum of the distributions are slightly different. This is due to the strong coupling of 
fo to kaons and, as already argued by Achasov [22j , the meson loop distributions coincide with the Flatte ones only 
in the limit of weak coupling. 

In Fig. [7] we show the decay rates fo — * 7T7t, fo — * KK and ao — > 7rry, ao — > KK as function of gj oK -^- The dashed 
areas in both plots correspond to the total decay rate of fo and ao as indicated by the PDG (notice, however, that in a 
note in PDG it is specified that the real width could be larger). To get agreement between our theoretical total decay 
rates and the measured ones, 9t Q j{K nas to lye between 3-4 GeV. The outcoming branching ratio ^ ao KK/^ a o^v ^ s 
- 0.3-0.4, larger than the PDG average 0.183± 0.024, while the obtained ratio T fo7T7r / (T fo7T7r + T. qK ^) is - 0.48-0.56 
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T[GeV] r[GeV] 




FIG. 7: Left panel: decay rates of fo as functions of the coupling constant g. The thin solid line corresponds to the decay into 
two pions, the dashed line to the decay into two kaons and the thick solid line is the sum of the two decay rates. 
Right panel: decay rates of ao as functions of the coupling constant g. The thin solid line corresponds to the decay into pion 
and eta, the dashed line to the decay into two kaons and the thick solid line is the sum of the two decay rates. 



in qualitative agreement with the (not bold) results listed in PDG. Notice furthermore that the fo mesons turns out 



to have typically a larger width than ao meson in agreement with Ref . 24 1 . 

We finally comment on a possible tetraquark unified interpretation of the light scalar mesons as presented in Ref. 
[Tl| . A too small decay constant g f oK j< would also imply a by far too narrow a and k mesons (related by Clebsh- 
Gordon coefficients [HI), thus against a tetraquark nonet. On the contrary, 9j oK -^ between 3-4 GeV is in agreement 
with a tetraquark nonet below 1 GeV, although problems, such as a too narrow k, persist, see discussions in [TlL fl5| . 
Such a strong coupling in the K K- channel implies that virtual cloud of kaon-antikaon pairs plays an important 

role, in particular for the fo resonance. A heuristic indicator of the mesonic cloud can be given by the quantity 

-l 



Z^ K = ( 1 + 2m \ dx K ) _ M J where R S ^ K = 9 2 s -j^ k R^kk(x)] (with S = fo,a ) refers to the kaonic loop 

only. As discussed in Section II in the subthreshold case (which applies to the kaonic channel here) the quantity 
(1 — Z) varies between and 1 and measures the mesonic cloud dressing of the original bare resonance S. In the fo 
case, by using Eq. (|5tj|) together with gf oK -^ = 3 GeV, one finds (1 — Z-g K ) — 0.38, hence implying a 38% of kaonic 
cloud. This number increases for increasing coupling strength gj gK j^- This discussion confirms the interpretation put 
forward in Ref. 0], where the light scalar mesons posses a tetraquark core but are dressed by kaonic clouds. 



IV. SUMMARY AND CONCLUSIONS 



In this work we studied the spectral functions of scalar mesons in one- and two-channel cases suitable for the de- 
scription of light scalar mesons below 1 GeV. We have computed, by using Lagrangians with non-derivative couplings, 
the propagators of scalar mesons at one-loop level. They satisfy for large ranges of parameters the Kallen-Lehmann 
representation, therefore implying normalized spectral functions. In this way a correct definition of decay amplitudes, 
weighted over the spectral function, is formulated: the finite-width effects are automatically taken into account. The 
resulting decay rates are smaller than the tree-level ones with increasing mismatch for increasing interaction strength. 
On the other hand, a sub-threshold tree-level forbidden decay, such as the KK mode for eto(980) and /o(980), becomes 
large. 

The resulting spectral functions for the a and k mesons may deviate consistently from the Breit-Wigner form. The 
Flatte distribution, although it approximates to a good level of accuracy the ao(980) and /o(980) spectral functions, 
emerges as a small-coupling limit of our more general spectral function. 

As stressed by Achasov [221 ] it is important to use distributions satisfying Kallen-Lehmann representations in 
experimental and theoretical studies. We thus believe that the use of distributions obtained from quantum field 
theoretical models fulfilling the correct normalization requirements can be helpful to correctly disentangle the nature of 
the scalar states. Future studies with derivative couplings, mixing effects and ^-decays represent a possible interesting 
outlook. 

Acknowledgements G.P. acknowledges financial support from INFN. 
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APPENDIX A: LOOP CONTRIBUTIONS 



Here we report basic formulas for the loop diagram of Eq. |(7j) drawn in Fig.l for the vertex function /a(?) = /a("9* 2 )- 
By evaluating the residua one obtains the one-dimensional integral 



s(x=p2) = i. /0 



dw- 



v w 2 + m 2 (4(w 2 + m 2 ) — x 2 ) 



(Al) 



which can be easily evaluated numerically for each well-behaved f\(w). We remind that within our conventions 

/a(0) = 1 and that w = k ; Eq. (|A1[) refers to a 3d-vertex function. In [22j the form f\(w) = 6 (A — w) is used and 

the limit A — > oo is taken. As described in the text, we did not follow this procedure but we used definite form(s) for 
the vertex function f\(w). As remarked in the text, we performed the calculations also with different forms for /a(ui) 
(different power form and exponential functions): the precise form of the cutoff function does not affect the physical 
picture. 

When the scalar state S couples to two particles of masses mi and mi the loop contribution is modified as following: 



w 



dw- 



yjw 2 + m\ + yjw 2 + m\ ^ f\ (w) 



\Jw 2 + m\ \Jw 2 + m\ ( {^\/w 2 + m\ + \/ w 2 + m 2 ^j 



(A2) 



The choice f\(w) = (l + w 2 /A 2 ) 1 with A = 1-2 GeV has been used in this work. 

In relation to the mass definition of Section II. B we report and compare in Table 1 the mass M defined in Eq. (I15p . 
the bare mass Mo, the maximum M max of the distribution ds=a(x) and the average mass (M) = J™dxxd s (x). We 
use m = m^, g a -K-K = 3 GeV, A = 1.5 GeV. 



Tab.l: Comparison of 'masses' 



M 


M 


M max 


(M) 


0.4 


0.61 


0.45 


0.54 


0.6 


0.71 


0.62 


0.65 


0.8 


0.86 


0.81 


0.82 


1 


1.03 


1.00 


1.01 



As expected, the larger the mass, the smaller the differences among the various mass-like quantities. 
APPENDIX B: SPECTRAL FUNCTION AS 'MASS DISTRIBUTION': AN INTUITIVE DISCUSSION 



We present an intuitive argument for the correctness of interpretation of the spectral function ds{x) as the 'mass 
distribution' of the state S. To this end we introduce two scalar fields A and B, the first massless and the second 
with Mb > Ms and write down the interaction Lagrangian 



L = cBAS + gSip 2 



(Bl) 



(for the following discussion the 'delocalization' of Eq. ((HJ is not important). We suppose that the interaction strength 
c is small enough to allow a tree-level analysis for the decay of the state B. The term cBAS generates the decay 
process B — ► AS, which reads (at tree-level) 



-vt-l 

BAS 



(M B ) 



8ttM 2[1 ' 



(B2) 



However, when g ^ the state S decays into (pip, that is the state S is not an asymptotic state. Physically, we observe 
a tree-body decay B — > Aipip, whose decay-rate reads: 



-it-l 

B Aipip 



(M B ) 



r 



t-i 

BAS 



(M B )d s (x)dx. 



(B3) 



The tree-body decay is decomposed into two steps: B — > AS and S — > ipip. The quantity ^bas^^b) represents the 
probability for B — > AS (at a given mass x for the state S) and ds(x)dx is the corresponding weight, i.e. the probability 



12 



that the resonance S has a mass between x and x + dx. In this example ds(x) emerges naturally as a mass distribution, 
correctly normalized, for the scalar state S. Furthermore, notice that in virtue of the limit d$(x) = 8{M — Ms) for 
g^O one has 

?ba w {Mb ) = r%\ s (M B ) for g - 0. (B4) 

In fact, if g is very small the state S is long-living and the equation (|B2[) is recovered. The present analysis also shows 
that studies on the tree-body decay of the meson can be consistent only if propagators satisfying the Kallen-Lehmann 
representation are used. 
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